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ON THE DESCENDING SERIES FOR BESSEL'S FUNCTIONS OF 

BOTH KINDS. 

By Pbof. James McMahon, Itb ' , N Y. 
1. It has been shown in the Annals* that Bessel's equation 

"^3+''^+('^-''')2/=o (1) 

has for complete primitive, when n is fractional, 

y = AJ„{x) + BJ_„ix) , (2) 

and when n is integral 

y = AJ„{x) + BK„ix) , (3) 

where 

fx^ 



JJx) 
and 



'2 



"""Y ( — 1)^ 



p=or(w + ;? + i)r(^ + i) 



ip 



(4) 



2^„((t') = ^^ + (-ir'^^. (5) 

Although this series for ^„(«), and that given on p. 88 for K,lx), are con- 
vergent for all values of x, yet when x is large the convergence is so slow as 
to render the series unavailable for numerical computation. 

In such a case the descending series are readily applicable ; and it is 
proposed to derive from equation (1) a general solution in negative powers of 
X, involving two arbitrary constants, and then to determine these constants for 
the cases of JJx), KJ^x), respectively. 

2. General solution in negative powers. 

In equation (1) put y = a^w, then 

a^^ + (2m. + l)x^+ O*^ + m^-n^)u = 0. 
To remove the second term, let m = — ^, then 

S + ^-e^^-i)^' («) 

in which u = x^y. If the right hand member were zero, the solution would be 

u = Ae*^ + Be-*" . (7) 

* Vol. VI, No. 4, p. 86. " On Bessel's Functions of the Second Kind," by Dr. Mazime Bocher. 
Annals of Mathematics, Vol. VIII, No. 2. 
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Substitute in (6), regarding A, B &8 functions of x, and identify separately the 
coefficients of «**, e~^ ; then 



(8) 



Put ix=^ 2; then 

^lA. 4- 2 ~ — Cn^ — it — 

d^ ^ ^ dz-^"" ^' ^■ 

Let A = A, + A,B-' + ... + A^-' + ..., (9) 

and substitute and equate coefficients of z~'~^ ; 

A _ A g(s + l)+i-n' _ . (2g + 1)" - 4n' 
»+'" ' 2(5 + 1) ~ ' 8(a-+1) ■ 

Substitute in (9) for A^, A^, . . . , and replace z by ix ; then 

A- Ah 4- 1'-^"' + (1' -4n')(y-4n') "I ] 

similarly, S- (10) 

.-. , from (7), 

x^y = A,e<^[l + K^ + ...]+ ^„a-[l - ^^Z.^.^ + ...], (11) 

in which A„, B„ are arbitrary constants. 

To make this expression real in form, put 

A, = iPe--', B, = iPe-', 
and then replace imaginary exponentials by trigonometric equivalents ; 

••• '"'y^ ^^'>^("' - "^L^- 1.2 (8^ — ^ ■ ■ J 

, „• / .rP — 4:n^ (P — 4?i2)(3^ — 4«2)(5='-4«2) , "l ,,0, 
+ Psm(«.-«)|_-^-g^-^ 1^;^ 3 (8^^ ^ + ...J.(12) 

This expression for y is a general solution of equation (1), involving the two 
arbitrary constants P, a. 
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3. To deterinine the constants P, a, in the case of Jjx). 
In the first place it may be verified that the ascending series for JJix) is 
equivalent to the definite integral 

a;"!" " 



= I I cos {x cos (p) H\v?^f df ; (13) 

-> 



l/;rr(«+i) 
for the coefficient of si?'' under the integral sign is 

t^^ Jcos^*^ 8in> df , 



which equals* 

{2p)\ r(n-i-p + i) • 

Then writing Pip + i) = (p - i) {p ~ i) ■ ■ ■ i PH) - -J^ l^i , and re- 
placing jo! hy r{p -f 1), the identity of the general term of (13) with that of 
(4) will be apparent. 

It is next proposed to determine the constants P, a, by comparing the 
forms taken by (12) and (13) when x is very large.t For convenience, let x be 
of the form 2r7i, where r is a large integer, then 

cos {x cos f) = cos (2a3 cos^ ^ (p — x) = cos (2a? cos^ ^ <p) , 
and the integral takes the form 

IT 

j cos (2a! cos^ ^ <p) sin^ip d<p , 



which, by the substitution 2x cos^ ^ f ^= t, becomes 

?^ r cos < . r-l (1 — i tx-')"-idt . 



X 



Taking this integral between and x, and then between x and 2a!, and trans- 
forming the latter part so that its limits become and x, it will be seen that 
the two parts are equal, and the integral becomes 



— ^, fcos t . r- i(l — \ tx-'^Y-^dt; 



* 'Williamson's Integral Calculus, Art. 122. 

t The form taken by (131 when x is large is suggested by Art. 406, Todhunter's Functions. 
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whence (13) ultimately takes the form 

211 



J^{x) = 






which reduces to 



\j:x 



cos i (w + i) K- .■* 



But the general solution for y in (12) becomes, on the same supposition (viz. 
that X is large and of the form 2rn-), 



P P 

y = —;= COS {x — a) = -— = cos a ; 
y X yx 



whence it follows that 



r2 
is 



COS a = cos i (m -\- ^)t: . 



(U) 



Again let x be of the form 2rn + ^n ; then cos (2a! cos^ \tp) wiU be replaced 
by sin (2a; cos^ ^(p), and the form of J„{x) will be 



00 

[;ra;j r(n + i) J 



which reduces to 
for y in (12) becomes 



' 2 1 1 . 
— I sin i (« + i) JT ; and, at the same time, the expression 



p p 

y = — = cos (x — «) = —7= sin a ; 

va; yx 



sin a ^ sin J (rt + ^) TT . 
From this and (14) it follows that 

whence 



Inx 
\2- 



_ J^{x) = cos (a; — J a- — ^ wh-) ^„(a!) + sin (a; — i «■ — ^ nn) (PJx)\ (15) 

for all values of x and n ; f„(a!), ^„(a!) stand for the descending series in (12). 

* Williamson's Int. Cal., Art. 124. 

t While this is a known result, I have not been able to find any previous complete derivation 
of it from Bessel's equation, with a general determination of the constants. 

An interesting application of the descending series for J^(x) is made by Sir G. G. Stokes, 
(Collected Papers, or Camb. Phil. Trans. Vol. IX) where he uses it to compute all the roots of the 
equation J^{x) = 0. Lord Bayleigh applies the same method to the equation Ji(x) = (Theory 
of Sound, Vol. 1, p. 273). 
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4. To determine the constants P, a, in the case of K^ix). 

Use (5), substitute for t/„(a;) from (15), and for convenience let x be any 
large number ; then 

J^ <4(«) = cos (je — i ff — -^ n^r) ; 

.-. 2 ^^ KJx) = I {sin (x-\n-\n7t)-^r (- 1)" sin (a; - i ;r + ^ nn)]. 

Write d for x — \-!z — \m: , then 

sin (a; — i '^ + i rm) = sin {d + nn) =^ ( — 1)" sin d , 
since, in this case, n is integral ; 

••• ^^n(-') = ^|sin^ = ^||cos(^-i.) 

= aIs^ ^^^ (x — ^Jz — ^nT:). 

\'2a! 

Comparing with (12), it follows that 

.-. J^ Klx) r- cos (« - f ;r - i nn) .pjix) + sin (« - J ;r - i m:) 4>„{x) * (1 6) 

for any value of x and an integral value of n ; (f„{x), </)„{x) stand, as before, for 
the descending series in (12). 

5. These series for J„{x), KJ^x) terminate only when 2n is an odd integer. 
In other cases their terms decrease in magnitude for a certain number of such 
terms (depending on the values of n and x), and then they begin to increase ; 
but the convergent part may be used for computation when x is large ; " for it 
can be proved that the sum of any number of terms of the series differs from 
the true value of the function by less than the last term included."t 

♦Equation (16) may be derived directly from (5) and (15) without previously assuming that 
Kn(x) must be capable of taking the general form given in (12). 

t Eayleigh, Theory of Sound, Vol. I, p. 264. Lommel, Bessel'sohe Punktionen, p. 62. 



